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The  article  presents  the  principles  of  investigation  of  physical 
fields  in  the  ocean  with  adaptation  of  the  measuring  algorithm  to 
the  statistical  structure  of  a field  and  describes  the  results  of 
experiments. 


An  important  condition  for  obtaining  reliable  scientific  data  on  the  physical 
processes  in  the  seas  and  oceans  is  the  matching  of  the  parameters  of  the  meters  and 
the  method  of  their  application  with  the  statistical  structure  of  the  hydrophysical 
fields  being  studied.  As  the  matching  criterion,  use  may  be  made  in  the  general 
case  of  the  relationship,  predetermined  with  a definite  degree  of  probability, 
between  the  frequency  band  of  the  meter  Af^  and  the  "spectral  window"  Af^,  in  which 

practically  the  entire  fluctuation  energy  of  the  investigated  element  of  the  hydro- 
physical field  is  concentrated. 


The  traditional  methods  of  research  in  the  sea  and  modern  technical  devices 
used  in  oceanography  do  not  at  the  present  time  permit  the  achievement  of  the  re- 
quired statistical  agreement  between  the  measurements  and  the  field  under  study. 

In  most  cases,  the  algorithms  of  measurements  conducted  in  the  sea  remain  unchanged 
during  the  period  of  observation  and  as  a rule,  the  instruments  operate  with  a 
fixed  speed  of  response  that  is  independent  of  the  structure  of  the  process.  In 
sounding  or  towing,  the  speed  of  the  relative  motion  of  the  sensor  and  medium  is 
also  maintained  constant  irrespective  of  the  dependence  on  the  configuration  and 
variability  of  the  field  element  being  measured.  This  method  of  investigation 
has  significant  drawbacks,  since,  as  is  well  known, distortions  of  measurement 
results  and  processing  of  primary  data  depend  not  only  on  the  characteristics  of 
the  meter  and  measuring  conditions  (dimensions  of  the  sensors,  their  orientation 
and  motion  in  the  medium  studied,  method  of  converting  the  signal  to  digital  form 
and  frequency  of  conversions,  duration  of  the  recording,  etc.),  which  during  the 
observation  remain  unchanged  and  are  known  in  advance,  but  also  on  the  structure 
of  the  investigated  process  itself.  Since  the  statistics  of  the  field  element 
studied  (and  hence,  Af  ) are  unknown  during  the  measurement,  and  a second  experi- 
ment under  identical  nsea  conditions  is  impossible,  there  is  always  the  following 
alternative: 


The  first  inequality  signifies  that  the  meter  gives  excess  data,  and  in  addi- 
tion, because  of  the  limitations  of  a series  of  observations  during  the  analysis, 
restricts  the  scale  to  the  range  of  low  frequencies.  Figure  1 a shows  a spectrum 
of  transparency  fluctuations,  the  frequency  of  whose  measurements  considerably  ex- 
ceeded Af  . The  "spectral  window"  was  "not  filled"  in  this  case.  The  fluctuation 
spectrum  nshows  a spike  in  the  range  of  lowest  frequencies  f0,  after  which  the 
spectral  intensity  retains  a low  constant  value  determined  by  the  meter  noise. 


Numbers  in  the  right  margin  indicate  pagination  in  the  original  text 


Figure  1.  a)  Spectrum  of  seawater  transparency  fluctuations,  bounded  In  the 
low-frequency  range;  b)  Spectrum  of  fluctuations  of  sea  surface  elevation, 
measured  with  large  distortions. 


When  Afy  < AfQ  (the  frequency  band  of  the  physical  process  is  greater  than 

that  of  the  meter),  the  results  of  the  study  are  distorted  by  the  suppression  or 
transformation  of  the  high-frequency  components  of  the  process.  An  example  of 
this  Is  the  fluctuation  spectrum  of  elevations  of  the  sea  surface  (Fig.  lb), 
obtained  by  analyzing  the  data  of  measurements  whose  frequency  was  smaller  than  fQ. 

In  the  same  plot,  the  dashed  curve  shows  the  dependence  — “ (y^ — )-^,  which 

corresponds  to  the  Phillips  spectrum  for  wind  waves.  As  Is  evident  from  the  figure 
the  fluctuation  spectrum  is  highly  distorted  and  does  not  permit  a correct  Inter- 
pretation of  the  result  of  the  measurements. 


Numerous  studies  have  shown  that  the  statistical  structure  of  most  real 
hydrophysical  fields  Is  variable  in  time  and  space.  The  unsteadiness  of  physical 
fields  In  the  ocean  Is  related  to  seasonal,  diurnal,  tidal,  and  other  oscillations 
In  some  cases,  an  intermittent  turbulence  is  observed,  when  periods  of  low  Inten- 
sity of  the  fluctuations  alternate  with  periods  of  their  appreciable  increase.^ 

In  all  these  cases,  the  value  of  Af  changes. 


The  lack  of  sufficiently  complete  a priori  Information  and  the  unsteadiness 
of  the  physical  process  do  not  permit  one  to  determine  in  advance  the  conditions  of 
measurement  and  characteristics  of  the  meters,  or  to  correct  the  results  of  the 
measurements.  Therefore,  In  order  to  limit  the  errors  of  measurement  and  processlni 
to  a predetermined  level  and  obtain  reliable  information.  It  Is  necessary  to  change 
the  algorithm  of  the  measurements  automatically  during  the  studies. 


Studies  characterized  by  a rearrangement  of  the  measurement  algorithm  and 
Its  "fitting"  to  the  statistical  structure  of  the  investigated  field  for  the  purpose 
of  obtaining  an  acceptable  and  constant  level  of  distortion  of  the  results  or  limiting 
the  lower  frequencies  are  called  adaptive,  and  the  measuring  devices  or  systems  which 
perform  this  rearrangement  In  the  real  time  of  the  measurements  are  called  adaptive 
meters  or  systems. 


The  fundamental  principles  of  operation  of  adaptive  meters  are:  analysis  of 
the  measured  signal  in  real  time  to  find  the  statistical  element  making  the  chief 
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contribution  to  the  error  of  measurement  and  processing;  comparison  of  this  para- 
meter with  a given  criterion  of  accuracy  and,  depending  on  the  result  of  the  com- 
parison, appropriate  modification  of  the  measurement  algorithm. 


Analysis  of  the  distortions  of  a random  field  by  real  meters  shows  ’ that 
the  chief  contribution  to  the  measurement  errors  is  made  by  small-scale  components 
corresponding  to  the  Nyqulst  frequency,  or,  in  other  words,  field  components 
whose  scale  is  comparable  to  or  smaller  than  the  minimum  scale  characteristic  of 
a specific  meter.  Therefore,  to  check  the  distortion  level  in  real  time,  it  is 
sufficient  to  have  information  on  the  behavior  of  the  spectrum  of  the  investigated 
process  in  the  high-frequency  region,  bounded  from  above  by  the  minimum  scale  of 
the  meter  Itself.  The  statistical  parameter  characterizing  the  behavior  of  the 
high-frequency  region  of  the  spectrum  is  determined  by  the  asymptotic  expansion  of 
the  Integral  Fourier  transform  of  the  correlated  function  R(t)  ■ a*p(r)  of  this 
process x 


where  s (w)  is  the  spectrum  of  the  process:  is  the  nth  derivative  of  the  corre- 
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lation  coefficient  when  t ■+•  0;  a is  the  variance. 


The  range  of  characteristic  dimensions  of  hydrophysical  meters  is  bounded 
from  below  by  values  on  the  order  of  fractions  of  a millimeter  and  hundredths  of 
a second,  and  from  above,  by  several  meters  and  seconds.  These  dimensions  correspond 
to  portions  of  the  spectrum  of  ocean  turbulence  located  in  the  region  of  universal 
equilibrium,  where  there  is  no  energy  input  to  the  fluctuations,  and  the  spectra 
are  characterized  by  a continuous  decrease  and  are  structurally  similar  to  Markov 
spectra. ^ A transformation  of  Eq.  (1)  for  such  processes  gives  the  following  ex- 
pression for  the  high-frequency  portion  of  the  spectrum  s (to) : 


where  D(A)  is  the  dispersion  of  increments  of  the  measured  process  during  the  minimum 
Interval  between  the  measurements  T,  . 


A formal  analysis  of  the  errors  of  reconstruction  of  the  initial  observation 
from  discrete  readings  and  of  the  errors  of  low-frequency  filtering,  calculation  of 
the  spectrum  and  conversion  of  the  signal  to  digital  form  by  various  digital-analog 
converters  shows  that  there  exists  a linear  or  quasl-linear  relationship  between 
the  dispersion  of  these  errors  and  D(A)  The  coefficients  of  proportionality 

between  them  depend  on  the  method  of  discretization,  modulus  of  the  frequency  char- 
acteristic of  the  sensor,  time  between  measurements,  and  other  parameters  known  in 
the  course  of  the  experiment. 


Thus,  the  dispersion  of  Increments  during  the  minimum  possible  Interval  between 
measurements  D(A)  is  a statistical  element  determining  the  behavior  of  the  high-fre- 
quency region  of  the  spectrum  and  hence,  the  measurement  and  data  processing  errors. 


The  dispersion  of  the  Increments  can  be  obtained  with  a given  probability  in 
the  real  time  of  observations,  this  being  a necessary  condition  for  fitting  the 
meter  to  the  structure  of  the  field  etudied. 


Analysis  shows  that  in  order  Co  determine  the  algorithm  of  adaptive  measure- 
ments, it  is  more  convenient  to  use,  not  the  sample  estimate  of  D(A),  but  the  mean 
modulus  of  Increments  M(|A|)  related  to  this  estimate.  This  relationship  for  non- 
Gauss  lan  random  processes  has  the  following  form: 

-fa)' 

Y, 

where  ■ y is  the  excess  coefficient  of  Increments  A of  a random  process  x(t) 

having  an  excess  coefficient  y . The  rate  of  convergence  of  the  sample  quantity 

x a 

|A|q  to  its  general  value  I ~~~  ) ' I dl{ * *(MI)  is  approximately  twice  as  high 

as  the  analogous  convergence  rate  of  the  sample  Increment  dispersion,  and  is  deter- 
mined from  a Cheyshev  inequality,  which  after  transformations  is  written  as 
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where  |A|q  Is  the  sample  value  of  M( | A ( ) at  n points;  1 - p is  the  estimation  proba- 
bility of  M(  | A | ) . Replacement  of  the  calculation  of  increment  dispersion  by  the 
mean  modulus  of  increments  in  real  time  appreciably  simplifies  the  technical  appli- 
cation of  the  adaptive  principle  of  measurement  and  makes  it  possible  to  reduce  the 
time  of  signal  analysis. 

It  follows  from  expression  (2)  that  for  a 20Z  spread  of  the  sample  values 
relative  to  the  true  value  (a  spread  of  this  magnitude  or  larger  takes  place  in 
the  study  of  the  statistical  structure  of  physical  fields  in  the  ocean*-*) , averaging 
can  be  performed  over  63  measurements  with  a probability  of  0.8.  The  notation  of 
the  algorithm  for  calculating  the  estimate  of  the  asymptotic  element  of  the  small- 
scale  fluctuation  component  model,  for  which  the  mean  modulus  of  increments  is  used, 
is  as  follows: 


m C ' 


where  m is  the  scale  factor  and  n is  the  averaging  step.  For  five  types  of  commonly 
used  digital  computers,  the  duration  of  the  computations  and  volume  of  the  storage 
device  Pv  were  determined  for  the  same  measuring  channel.  The  results  of  the 
calculations  are  given  in  the  table. 


Table 


Parameter 


Minsk-22 


Duration  of  calculations. 

Msec 

Volume  of  SD  for  calcula- 
tions, cells 

Total  volume  of  SD 


Type  of  computer 
Dnepr-1 


M-220 


BESM-4 


1460 

666 

520 

802 

68 

68 

68 

68 

8192 

1536 

4096 

4096 

If  the  number  of  measuring  channels  is  nv,  and  the  maximum  frequency  of 
measurements  in  each  channel  is  f , the  boundary  relations  for  the  computer  caps' 

v i 

billtles  may  be  written  as  follows : — > n f , P * n P , where  P is  the 

v v comp  v v comp 

volume  of  the  storage  cell  of  a given  computer  type. 


Figure  2.  Region  of  possible  use  of  different  types  of  digital  computers  for  cal' 
culatlng  in  real  time  the  asymptotic  element  of  the  model  of  the  small-scale  com- 
ponent of  fluctuation  processes. 


On  the  basis  of  the  above  relations,  Fig.  2 shows  the  regions  of  possible 
application  of  different  computer  types  as  calculators  of  a multichannel  adaptive 
measuring  system.  The  plots  show  that  for  a complete  loading  of  the  computer  at 
a maximum  measuring  frequency  of  20  times/sec  in  each  channel,  the  limiting  number 
of  channels  for  different  computers  ranges  from  22  to  32.  This  indicates  that 
universal-type  computers  should  be  used  only  for  adaptive  control  of  meters  inte- 
grated in  a single  automated  system  and  operating  simultaneously.  For  individual 
meters  in  which  the  number  of  channels  is  no  greater  than  10,  these  types  of  cal- 
culators will  not  be  used  effectively.  In  this  case,  it  is  preferable  to  entrust 
the  calculation  of  the  asymptotic  element  of  the  small-scale  component  model  of  a 
hydrophysical  process  to  a specialised  calculator  connected  to  the  feedback  circuit 
of  the  measuring  unit  of  the  equipment. 


On  the  basis  of  the  above  principles,  the  Marine  Hydrophysical  Institute  of 
the  Ukrainian  Academy  of  Sciences  has  designed  a measuring  unit  of  the  equipment 
for  adaptive  studies  of  physical  processes  in  the  sea.  This  unit  was  tested 
on  an  experimental  sea  range.  The  general  block  diagram  of  the  set  is  given  in 
Fig.  3.  The  equipment  includes:  a set  of  sensors  of  hydrophysical  characteristics 
(D^,  D2,...,Dn),  a multichannel  dc  amplifier  unit  (A),  and  device  for  receiving, 

controlling  the  recording  and  coding  information  (CU),  a 15-track  tape  recorder  (TR), 
a set  of  electronic  analog  recorders  (ER.  - ER.) , a digital  printer  (DG),  an  electric 
puncher  (EP),  and  an  adaptive  control  device  iAD). 


i 
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Figure  3.  Block  diagram  of  measuring  unit  for  adaptive  measurements:  a)  under 
measuring  conditions;  b)  under  reproduction  conditions. 


The  newly  designed  equipment  was  used  to  measure  the  elevation  of  the  sea 
surface,  wind  pressure  on  the  wave,  horizontal  components  of  the  flow  velocity  on 
several  levels,  and  fluctuations  of  the  index  of  attenuated  directed  light  and 
temperature.  The  unit  was  designed  for  synchronous  operation  on  8 analog  and  8 
digital  measuring  channels,  with  a measuring  frequency  of  25  Hz  in  each  channel. 

The  recording  equipment  can  operate  in  the  cyclic  or  slave  mode.  The  hydrophysical 
parameter  sensors  were  Installed  on  a stationary  gradient  mast  rigidly  mounted  on 
a concrete  support  at  a depth  of  15  m and  a distance  of  300  m from  the  shore. 

The  measuring  signals  from  the  sensor  outputs  were  transmitted  to  a shore  labora-  /102 
tory  over  a multicore  cable, and  after  passing  through  a multichannel  amplifier, 
reached  the  Inputs  of  a control  unit  (CU)  and  an  adaptive  control  device  (AD). 

The  latter'  Included  an  analog-code  converter,  a memory  register,  a computer,  and 
a system  for  automatic  control  of  the  measurements. 


The  output  signal  of  the  amplifier  is  converted  to  the  digital  form  Y(1T.  ) 
with  period  T^  between  readings.  The  digital  data  are  fed  to  the  computer,  where 
the  sample  value  of  the  mean  modulus  of  Increments  between  neighboring  numbers  is 
determined  In  real  time  from  the  given  averaging  step.  The  result  of  the  calcula- 


tions of  |a|q  ■ — £ [Y^  - y1+1],  whose  magnitude  is  determined  by  the  structure 

of  the  physical  process  being  measured,  enters  the  input  of  an  automatic  control 
system,  where  it  is  compared  with  a preassigned  value  of  the  mean  modulus  of  incre- 
ments M(Tfc),  which  is  related  to  the  allowed  distortion  level  of  the  auto-  and  mutual 
spectra  characterizing  the  process  studied.  Depending  on  the  sign  and  magnitude  of 

the  difference  of  J^T^)  and  |a|q,  a signal  acting  on  the  measurement  algorithm  Is 

formed  at  the  output  of  the  system.  In  the  elaborsted  unit, this  signal  controls 
the  oscillator  frequency  synchronizing  the  operation  of  the  analog-code  converter 
and  channel  commutator  of  the  control  unit,  and  thereby  performs  the  sdaptlve  matching 
of  the  pas abend  of  the  meter  with  the  "spectral  window"  of  the  process  studied. 


The  value  of  Che  measurement  frequency  is  recorded  In  analog  form  with  an 
electronic  plotter  and  In  digital  form  with  an  electronic  counting  frequency  meter 
The  output  voltage  from  the  comparison  unit  of  the  control  system  may  be  used  to 
control  the  measurements  In  other  channels,  for  example,  the  rate  of  Immersion  of 
a towed  or  sounding  device,  orientation  of  the  sensor  In  the  flow,  dimensions  of 
the  measuring  base  of  the  Instrument  or  system  of  sensors,  etc. 


The  basic  condition  for  carrying  out  studies  by  means  of  an  adaptive  measuring 
unit  was  to  maintain  at  a given  level  the  relative  magnitude  of  distortions  (q)  of 
the  power  spectrum  of  fluctuation  processes  in  the  surface  layer  of  the  sea.  Irrespec- 
tive of  the  frequency  distribution  of  energy 

amax 


where  s is  the  maximum  value  of  the  spectrum  and  P_.  is  the  absolute  distortion 

level.  n 


The  quantity  q determines  the  decrease  of  the  spectrum  into  the  region  of 
small-scale  processes  and  may  be  given  by  taking  into  account  the  minimum  scale 
of  the  fluctuations  accessible  for  study.  The  scale  is  determined  by  the  struc- 
tural and  metrological  characteristics  of  the  meters. 


The  analysis  performed  showed  that  the  reciprocal  of  the  relative  error  of 
spectral  distortions  is  related  to  the  dispersion  of  increments  during  the  minimum 
time  between  measurements,  to  the  power  exponent  of  the  spectral  decrease,  to  the 
increase  in  frequency  and  to  the  dispersion  of  fluctuations  or  their  dynamic  range 
A , and  is  expressed  by  the  following  relation: 


where  is  the  variance  of  fluctuations  of  the  process  studied;  n , s are  the  power 

exponents  of  the  increasing  and  decreasing  branches  of  the  spectrum,  respectively; 
D(A)  is  the  dispersion  of  the  signal  increments. 


Since  the  end  result  of  the  studies  was  the  determination  of  the  spectral 
characteristics  of  the  processes,  it  was  convenient  to  assign  the  distortion 
level  in  terms  of  the  number  of  orders  P by  which  the  spectrum  decreases  from  the 
maximum  to  the  minimum  value.  From  the  ’latter  expression  it  is  easy  to  obtain 
Pa  - -log  q. 


where  A - 1.95-2 


Since  the  dispersion  of  increments  D(A)  is  uniquely  related  to  the  mean  modulus 
of  increments  M(|A|),  and  on  the  basis  of  the  Chebyshev  inequality  the  fluctuation 
dispersion  is  expressed  with  reliable  probability  (1  - P)  in  terms  of  the  magnitude 

y 2<JX 

of  the  dynamic  range  JL.  < , the  above  relation  may  be  written  as  follows: 


For  physical  processes  in  ocean-atmosphere  boundary  layers,  the  exponent  s 
of  the  spectral  decrease  lies  within  the  range  of  4-5,  so  that  the  preceding  ex- 
pression for  (1  - P)  ■ 0.8  may  be  reduced  to  the  following  simple  form: 


For  physical  processes  structurally  similar  to  a locally  isotropic  turbulence 
, the  relation  for  the  mean  modulus  of  increments  is 


For  a given  Pq  and  hence  MCT^) , the  adaptive  meter  calculates  in  real  time 

the  mean  modulus  of  increments  | A | , compares  it  with  M(T^) , and  automatically  finds 
the  measurement  frequency  required,  maintaining  it  at  a constant  level,  if  the 
measured  physical  process  has  properties  similar  to  those  of  a stationary  one. 
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Figure  4.  Typical  trajectory  of  adaptive  scanning  of  measurement  frequency 


Figure  4 shows  a typical  scanning  trace  of  the  measurements  with  an 
adaptive  meter.  Initially,  when  the  feedback  circuit  was  disconnected , the  measure- 
ment frequency  was  maximum  and  equal  to  15  Hz.  At  the  Instant  when  the  feedback 
was  connected,  the  measurement  frequency  decreased  to  a certain  level  that  was  sub- 
sequently kept  constant  with  only  slight  fluctuations  if  the  process  was  stationary 


Figure  5 shows  elevation  spectra  of  the  sea  surface  obtained  with  an  adaptive 
meter  under  various  hydrometeorological  conditions  and  a given  value  P-  - 3.  In  this 

f "1 

figure,  the  spectra  were  plotted  in  normalized  form:  1 - 2m  — is  laid  off  along  the 


abscisaa  axis,  and  , along  the  ordinate  axis.  The  graph  shows  that  a specified 

max 

spectral  decrease  of  three  orders  of  magnitude  is  provided  by  the  meter.  The  mea- 
surement frequency,  whose  magnitude  was  established  by  the  adaptive  meter  in  each 
specific  case,  changed  appreciably,  as  is  evident  from  the  figure.  A developed 
wind  wave  at  a wind  speed  of  7.5— 8.5  m/sec  and  P.  * 3 is  associated  with  a measure- 
ment frequency  of  9-10  Hz;  calm  weather  with  a wind  speed  of  1. 0-2.0  m/sec  and  P - 3 
is  associated  with  a frequency  of  1.5-2  Hz.  q 

Simultaneously  with  adaptive  measurements,  traditional  measurements  were  also  /105 
performed,  in  which  the  observation  frequency  was  4 Hz  regardless  of  the  structure 
of  the  process  studied.  The  results  of  these  measurements  (Fig.  6)  show  that  the 
relative  distortion  level  may  vary  by  a factor  of  10  under  different  conditions. 

Adaptive  methods  with  a changing  measurement  frequency  make  is  possible  to  /106 
keep  constant  the  smallest  investigated  scale  of  the  process  in  unsteady  changes 
of  its  structure.  This  fact  becomes  particularly  important  in  the  study  of  the  /107 

statistical  properties  of  ocean  fields  with  well-defined  unsteady  properties  such 
as  intermittent  turbulence,  irregularity  of  the  intensity  of  turbulent  pulsations 
with  depth  in  a stratified  ocean,  etc. 
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